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A  self  contained  program  has  been  developed  to  meet  the  objec¬ 
tive:  (1)  Generation  of  thermochemical  properties  for  any  assigned 
propellant ,  (2)  Analysis  of  transient,  inviscid,  compressible  flow 
through  the  gun  barrel,  (3)  Analysis  of  unsteady,  viscous,  compressible 
flow  with  arbitrary  pressure  gradients  on  the  bore  surface,  (4)  Analy¬ 
sis  of  transient  heat  diffusion  through  single  or  multi-layer  gun  tube, 
(5)  Analysis  of  unsteady  free  convection  and  radiation  around  gun  uube, 
and  (6)  Control  of  sequence  of  procedures  due  to  the  nature  of  unknown 
and  transient  boundary  conditions  and  also  due  to  coupH«g  between 
problems.  The  investigations  in  these  areas  revealed:  Carbon  monoxide 
is  more  than  40%,  supersonic  velocity  it  the.  base  of  the  projectile 
for  a  constant  diameter  tube  even  though  the  gas  temperature  is  high 
and  also  Impossible  for  steady  state  flows, period  in  which  LaCranglan 
assumption  may  be  approximate  enough,  extension  of  the  state  of  the 
art  in  unsteady  viscous  flows,  need  for  new  approaches  without  the  use 
of  analogy  between  momentum  and  energy  transfer,  need  for  inclusion  of 
non-Fourier  effects  in  classical  Fourier  heat  conduction  models,  and 
order  of  magnitude  for  free  convection  and  radiation.  It  is  obvious 
that  there  is  strong  coupling  between  unsteady  viscous  flow  and  un¬ 
steady  heat  diffusion.  For  hypervelocity  systems,  there  is  a  possi¬ 
bility  of  interaction  between  inviscid  and  viscous  flows  as  in  hyper¬ 
sonic  flow.  When  all  these  programs  are  arranged  for  a  digital  com¬ 
puter,  the  overall  program  can  be  big  like  NASA's  NASTRAN  program. 

When  these  two  programs  are  utilized,  then  only  the  resulting  gun 
would  be  realistic. 
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INTRODUCTION 

Rapid  fire  automatic  weapons  permit  discharge  of  hundreds  of 
projectiles  per  minute.  As  a  result,  the  surface  temperatures  of  the 
gun  tube  easily  reach  1500°F  at  the  inner  surface,  and  1000°F  a _  the 
outer  surface  in  a  short  period  of  time,  even  though  both  surfaces  were 
initially  at  room  temperature.  There  is  only  one  approach  apparent 
that  can  fulfill  the  objective  stated  in  the  title  other  than  cut  and 
try  techniques.  Any  experimental  approach  has  limited  value  due  to 
moving  action  of  the  projectile  and  the  presence  of  large  transient 
gradients  at  the  bore  surface.  If  the  instantaneous  effects  are  not 
considered,  structural  failure  may  result  because  of  excessive  thermal 
stresses  near  the  surface  of  the  solid  where  structural  strength  was 
already  reduced  by  the  high  temperature. 

What  is  the  state  of  the  art  in  the  development  of  weapon  or 
weapon  systems?  It  is,  in  the  author’s  opinion,  a  cut  and  try  pro¬ 
cedure.  That  is,  in  other  words,  a  capability  does  not  exist  to  gen¬ 
erate  the  pertinent  information  regarding  the  behavior  of  the  weapon 
subjected  to  any  ammunition  and  firing  schedule  without  preliminary 
design,  prototype  development,  tests,  and  analyses.  The  next  logical 
question  is:  Is  it  possible  to  formulate  a  program  to  obtain  future 
weapon  characteristics  before  actual  development?  The  answer  may  be 
stated  as  "yes".  NASA’s  NASTRAN  program  is  quite  familiar  to  the 
structural  analysis  community  and  also  to  the  designers.  If  a  program 
can  be  developed  to  deal  with  the  influences  of  fluid  dynamic  aspects 
of  weapon  systems,  this  program  in  addition  to  NASTRAN  can  make  it 
possible  to  develop  future  weapon  systems  without  cut  and  try  tech¬ 
niques.  So,  the  rest  of  this  paper  will  deal  with  the  propellant 
gases.  In  particular,  its  influence  on  unsteady  heat  transfer  analysis 
for  any  ammunition,  gun,  and  firing  schedule. 

As  the  projectile  moves  ahead  because  of  the  high  pressure 
gases  created  by  burning  propellant,  the  propellant  gas  will  be  set 
into  motion  starting  from  rest.  Since  the  governing  equations  of 
fluid  dynamics  for  many  problems  of  interest  are  a  system  of  nonlinear 
partial  differential  aquations  and  are  also  dominated  by  real  gas  and 
nonequilibrium  effects,  no  general  solutions  exist  that  allow  arbitrary 
initial  and  boundary  conditions.  Therefore,  examination  of  the  flow 
field  and  subdivision  of  the  overall  problem  by  consideration  of  domi¬ 
nant  features  only  seems  appropriate.  The  objective  is  to  establish  a 
capability  to  perform  overall  heat  transfer  analysis  for  any  given 
dimensions  of  the\weapon  and  for  specified  propellant  characteristics. 
Towards  this  goal ,vthe  propellant  gas  convective  heat  transfer  problem 
is  divided  into  several  subproblems  f-l-J-j  2  -  Generation  of  thermo¬ 
chemical  properties  for  any  chosen  propellant;  11-*' Analysis  of 
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transient,  inviscid,  compressible  flow  through  the  gun  barrel;  HI  - 
Analysis  of  unsteady,  viscous,  compressible  flow  with  arbitrary  pres¬ 
sure  gradients  on  the  bore  surface;  TV--  Analysis  of  transient  heat 
diffusion  through  single  or  multi-layer  gun  tubes;  V--»  Analyst  of 
unsteady  free  convection  and  radiation  around  the  gun  tube;  anef-4#£<*“ 
Control  of  a  sequence  of  procedures  due  to  the  nature  of  unknown  and 
transient  boundary  conditions  and  also  due  to  coupling  between  the  pro¬ 
blems.  ,  Let  us  review  each  and  every  one  of  these  problems  in  the 
following  sections. 

I  THERMOCHEMISTRY 

The  generation  of  thermochemical  properties  for  any  propellant 
under  consideration  is  essential  to  the  development  of  a  capability  to 
perform  overall  heat  transfer  analysis  and  thereafter  stress  analysis 
for  any  given  dimensions  of  the  weapon  and  for  any  specified  propellant 
characteristics.  The  thermochemistry  of  propellants  involves  deter¬ 
mination  of  chemical  composition  of  propellant  gases  either  by  finite- 
rate  chemistry  or  by  chemical  equilibrium  chemistry  and  the  derivation 
of  propellant  gas  properties  from  the  composition  [2],  The  formulation 
of  a  finite-rate  chemistry  model  involves  the  following  events  [31: 

*  Selection  of  only  important  species,  such  as  a  dozen.  The  proper 
selection  may  depend  upon  the  actual  experiments. 

•  Selection  of  only  important  reactions.  Of  course,  the  number  of 
species  have  to  be  matched  with  the  number  of  reactions.  Again,  the 
proper  selection  may  depend  upon  the  experience,  intuition,  and  ob¬ 
servations. 

•  Determination  of  forward  and  backward  reaction  rates  -  sometimes, 
this  amounts  to  either  conducting  special  experiments  or  estimating 
them  from  kinetic  theory  of  gases  and  empirical  models. 

•  Formation  of  nonlinear  ordinary  differential  equations  for  rate 
of  change  of  species  -  one  for  each  specie. 

*  Linearization  (implies  no  significant  change  of  species  concen¬ 
trations  within  a  time  step)  by  Taylor's  series  expansion. 

*  Development  of  special  solution  procedures  to  solve  a  large  system 
of  linearized,  nonhomogeneous ,  first  order  ordinary  differential  equa¬ 
tions,  such  as,  combination  of  subdomain  method,  i.e. »  particular  case 
of  method  of  weighted  residuals  and  matrix  solution  procedure. 
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•  Postulation  of  "quasi-global"  mechanism  to  bridge  the  gap  between 
the  original  fuel  and  the  equilibrium  products  of  combustion  -  This 
step  is  essential  because  the  detailed  reaction  mechanism  becomes  ex¬ 
tremely  complex  even  for  hydrocarbon  fuels,  such  as,  kerosene  (C9H20). 
Of  course,  this  must  be  compatible  with  burning  rate  and  ignition  delay 
period. 


It  is  clear  that  the  above  events  up  to  a  certain  extent  are 
hypothetical  and  also  require  prior  knowledge.  The  feasibility  of 
such  a  formulation  is  definitely  questionable  for  any  new  (would  be) 
propellants.  It  is  important  to  look  into  the  structure  of  propellants 
and  the  available  information  on  such  propellants.  The  constituents 
of  typical  propellants  are  as  follows: 


Physical 

Chemical 

IMR 

M18 

Nitrocellulose 

C6  H7. 36S  °10.27  N2,635 

.8984 

.7960 

Nitroglycerin 

C3  H5  Og  Nj 

- 

.0995 

Dlbutylphthalate 

Cl6  H22  O4 

- 

.0895 

Dinitrotoluene 

c7  h6  o4  n2 

.0718 

- 

Diphenylamine 

Cj2  Hn  N 

.0062 

.0099 

Ethyl  Alcohol  (Residual) 

c2  h6  0 

.0053 

.0049 

Potassium  Sulfate 

k2  S04 

.0089 

- 

Water  (Residual) 

H20 

.0089 

- 

Since  the  Important  constituents  are  quite  complex  in  their 
structure  and  the  chain  of  reactions  leading  to  lower  compounds  is  un¬ 
known,  it  is  impossible  to  consider  finite-rate  chemistry  at  the  pre¬ 
sent  time.  However,  the  chemical  equilibrium  chemistry  can  be  con¬ 
sidered  with  possible  hundreds  of  species  and  reactions  in  ovder  to 
obtain  the  thermochemlcal  properties  of  typical  propellant  gases. 

The  information  on  chemical  equilibrium  compositions  of  a 
chemical  system  permits  one  to  calculate  theoretical  thermodynamic  and 
transport  properties  for  the  system.  Therefore,  the  first  step  is  to 
find  chemical  equilibrium  composition  for  any  assigned  thermodynamic 
state.  Chemical  equilibrium  is  usually  described  by  either  of  two 
equivalent  formulations,  l.e.,  equilibrium  constants  or  minimization 
of  free  energy.  However,  the  minimization  of  free  energy  is  somewhat 
better  and,  therefore,  it  is  utilized. 

The  condition  for  equilibrium  may  be  stated  in  terms  of  any  of 
several  thermodynamic  functions  such  as  the  minimization  of  the  Gibbs 
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free  energy  or  Helmholtz  free  energy  or  the  maximization  of  entropy. 

If  one  uses  temperature  and  pressure  to  characterize  a  thermodynamic 
state,  the  Gibbs  free  energy  is  most  easily  minimized  in  as  much  as 
temperature  and  pressure  are  its  natural  variables.  Similarly,  the 
Helmholtz  free  energy  is  most  easily  minimized  if  the  thermodynamic 
state  is  characterized  by  temperature  and  volume  or  density.  The  com¬ 
puter  program  based  on  these  formulations  is  utilized  to  predict  the 
composition  for  any  propellant  and  for  any  assigned  thermodynamic  state. 

The  theory  for  deriving  transport  properties  of  gases  from  in- 
termolecular  potentials  which  determine  the  forces  between  any  two 
molecules  in  the  gas  is  available  [4].  However,  the  determination  of 
intermolecular  potentials  between  complex  molecules  is  almost  im¬ 
possible.  Measurement  and  tabulation  of  transport  properties  such  as 
thermal  conductivity,  X,  and  dynamic  viscosity,  y,  for  all  gas  mixtures 
and  under  all  circumstances  is  not  feasible.  Thus,  a  means  of  pre¬ 
dicting  the  transport  properties  from  known  data  is  necessary.  Since 
computed  composition  of  gases  may  not  be  that  accurate  either,  the 
transport  properties  are  computed  by  summing  the  products  of  individual 
mole  fractions  and  corresponding  transport  properties.  The  individual 
transport  properties  at  the  desired  tremperature,  T,  are  obtained  by 
the  use  of  modified  Sutherland's  relation. 

Some  of  the  highlights  of  typical  propellant  gases  are  as 

follows: 

•  About  73  products  of  combustion  are  posaible  for  IMR.  However, 
the  Important  species  are:  CO  »  0.43,  fi2  *  0.12,  H20  -  0.21,  N2  ■  0.12, 
C0a  -  0.12 

•  About  48  (less  than  73,  due  to  absence  of  potassium  and  sulphur 
elements)  products  of  combustion  are  possible  for  M18.  The  important 
species  are:  CO  »  0.41,  H2  -  0.19,  H20  »  0.16,  H2  -  0.1,  C02  ■  0.08, 

•  The  gases  were  found  to  be  highly  toxic.  The  resultant  gases  may 
be  due  to  the  availability  of  lass  oxidants  in  those  propellants.  The 
consequences  of  Incomplete  combustion  are  nuzzle  flash,  smoke,  and  fire 
in  addition  to  low  efficiency.  '!oyever,  there  may  be  other  effects  that 
may  be  favorable  from  another  vw«  rede  ‘ .  The  ideal  propellant  gas 
should  possess  low  burning,  high  impetus,  and  low  molecular  weight 
characteristics  from  the  thermodynamic  point  of  view.  However,  the 
Impetus  is  directly  proportional  to  the  flame  temperature  end  Inversely 
proportional  to  the  molecular  weight.  Also,  the  elements  nitrogen  end 
oxygen  are  heavier  than  carbon  and  hydrogan.  Therefore,  the  Ideal 
thermodynamic  propertlae  may  ba  obtained  at  the  expense  of  laaa  end 
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less  oxidants  In  the  formation  of  propellants  and  also  by  promoting 
incomplete  combustion. 

•  The  major  component  of  present  day  propellants  is  nitrocellulose, 
c6H7.36sOio.27N2. 365  (molecular  weight  *  274.575,  oxygen  »  59.85  per 
cent).  Such  propellants  are  expected  to  yield  about  42  per  cent  of 
carbon  monoxide.  With  the  use  of  other  major  elements,  such  as,  nitra- 
mine-HMX  (C4H8N808)  being  investigated  by  other  elements  of  DOD,  such 
as  the  Air  Force,  the  carbon  monoxide  is  expected  to  be  more  than  50 
per  cent  of  propellant  gases.  The  oxygen  content  of  this  propellant 
(molecular  weight,  296)  is  only  43.24  per  cent.  However,  this  type  of 
propellant  showed  higher  burning  rates  and  a  larger  spike  in  the 
pressure-time  trace  (a  safety  hazard),  and  thus  contributed  to  lower 
ballistic  efficiency.  If  the  burning  rates  can  be  controlled  to  yield 
approximately  constant  pressures  in  a  chamber,  many  of  the  weapon  pro¬ 
blems  that  exist  today  can  be  solved. 

•  The  thermodynamic  properties  including  adiabatic  flame  temperature 
and  composition  are  determined  for  various  pressures.  The  variations 
with  pressure  is  Interpreted  as  due  to  variations  in  composition  with 
pressure.  The  Increase  in  pressure  causes  increase  in  molecular  weight, 
adiabatic  flame  temperature  and  specific  heats,  but  a  net  decrease  in 
specific  heats  ratio. 

•  The  transport  properties  of  propellant  gases  are  much  more  impor¬ 
tant  than  the  thermodynamic  properties  as  far  as  forced  convection  is 
concerned.  The  convective  heat  transfer  coefficient  is  directly  pro¬ 
portional  to  the  thermal  conductivity  of  gaaes.  Since  the  mean  free 
path  of  gaa  molecules  i©  inversely  proportional  to  the  pressure,  and 
number  of  molecules  per  unit  volume  is  directly  proportional  to  the 
pressure,  the  thermal  conductivity  (of  gasea)  which  ie  a  product  of  the 
two  Is  not  a  strong  function  of  pressure.  However,  the  thermal  con¬ 
ductivity  depends  strongly  on  the  temperature  of  gases.  The  thermal 
conductivity  increases  slightly  faster  than  the  dynamic  viscosity,  but 
increases  many  fold  (more  than  three  times,  500  *  2000°K  range)  with 
increase  in  temperature.  However,  the  Prandtl  number  ie  nearly  inde¬ 
pendent  of  temperature.  It  is  to  be  noted  that  the  variation  of  speci¬ 
fic  heat  with  temperature  ie  included  In  the  calculation  of  the  Prandtl 
number. 


II  CORE  FLOW 

Since  the  composition  and  thermodynamic  properties  are  esta¬ 
blished  in  the  above  section  for  any  chosen  propellant,  it  is  now 
logical  to  define  the  transient  inviecid  compressible  flow  through  the 
central  section  of  the  barrel.  Almost  all  predictions  of  the  behavior 
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of  propellant  systems  are  based  on  an  admixture  of  intuition,  experi¬ 
ence  and  the  use  of  interior  ballistic  equations.  These  equations 
assume  space-mean  averages  for  the  thermodynamic  variables  and  include 
dynamic  effects  only  through  the  inclusion  of  various  correction  fac¬ 
tors  determined  more  or  less,  empirically.  It  is  also  common  to  assume 
either  uniform  density  or  a  linear  velocity  gradient  in  the  flow  field 
between  the  breech  and  the  bullet  at  any  time  during  weapon  firing. 


The  core  flow  model,  while  necessarily  a  gross  over-simplifica¬ 
tion  of  the  physical  system,  would  permit  at  least  a  study  of  the  distri¬ 
bution  of  thermodynamic  variables  and  provide  more  accurate  Interior 
ballistics  data  than  the  data  obtained  by  conventional  methods  to  aid 
in  the  design  of  muzzle  brakes,  noise  and  flash  suppressors,  and  also 
gas  ports  for  automatic  actuating  mechanisms. 

A  model  involving  unsteady,  one-dimensional  motion  would 
apparently  be  the  most  that  could  be  tolerated  in  complexity.  More¬ 
over,  this  model  does  not  warrant  any  more  complexity  due  to  the 
limited  knowledge  in  solid  propellant  particle  aovements  and  in  asso¬ 
ciated  burning  models.  The  gas  is  assumed  to  be  inert  and  also  with¬ 
out  viscosity  and  heat  conduction  effects.  This  is  in  compliance  with 
Frandtl's  idea  of  separating  the  flow  into  core  flow  and  boundary 
layer  flow.  For  convenience,  the  mass  and  energy  are  assumed  to  be 
distributed  continuously  throughout  the  gas  flow,  although  such  an 
assumption  is  unnecessary  to  apply  the  chosen  method  of  approach  for 
the  solution  to  the  problem.  The  governing  equations  are  as  follows: 


Continuity: 


Momentum: 

3u  ..  ..  3u  k  "  1  3i 
P 


4*  U  1C~‘  *******  «•  Q 

at  ax  p  3x  p H 


D 


Energy: 


P  ^  («  +  f  )  +  P«  (•  ♦  §  )  +  Q  (•  +  -  )  <l*u> 

State: 
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The  wall  friction  does  not  appear  directly  in  the  energy  equa¬ 
tion  since  its  action  is  that  of  converting  mean  kinetic  energy  to 
thermal  energy.  The  above  governing  equations  of  unsteady,  inviscid, 
and  compressible  flow  are  hyperbolic  in  nature.  In  principle,  one  can 
solve  them  by  the  well-known  method  of  characteristics.  Along  certain 
(characteristic)  directions,  the  coupled  partial  differential  equations 
can  be  reduced  to  a  set  of  simultaneous  ordinary  differential  (charac¬ 
teristic)  equations.  The  solution  of  such  a  set  of  equations  is  fur¬ 
ther  complicated  by  an  unknown  nonhomogeneous,  and  unsteady  (moving 
bullet)  boundary  condition  and  time-dependent  propellant  burning  model. 


The  governing  equations  can  be  conveniently  written  with  u,  a 
(speed  of  sound),  and  s  (entropy)  as  dependent  variables  Instead  of  u, 
p  (density),  and  e  (internal  energy).  If  the  definition  of  perfect 
differentials  is  Included  as  additional  equations,  one  can  obtain  any 
partial  derivative  of  these  dependent  variables  by  Cramer's  rule  of 
determinants.  Since  the  characteristic  directions  can  be  defined  as 
the  curves  along  which  these  derivatives  are  discontinuous,  one  can 
obtain  the  following  sets  of  equations  by  forcing  respectively  the  de¬ 
nominator  and  the  numerator  of  any  one  of  those  partial  derivatives  to 
zero. 


Characteristic  directions: 


dx 

dt 


Du 


•  x  +  a 


dx 

dt 


•  x  -  a 


dx 

dt 


Governing  equations  along  characteristic  directions: 


Dt 


likm.  ■  JL  +  aQ(l-pn)(W  +  v*  -  Ct>T)/Cn 

Y-l+2pn  Dt  Cp  Y-1+2PH  Dt  2  ^ 

+  eQU/p-1/6)  -  uQ/p  -  2ful/D 

-  ACIlSn)..  Si  m  «  JL  iVrMUAnl,  Si  -  aQ(l-0r»)  (W  +  -  CnTVCn 

Dt  Y-l+lpn  Dt  Cp  Y**l+2pn  Dt  2 

-  aQ(l/p~l/6)  -  uQ/p  -  2fua/D 


Da 

Dt 


fit  [B  +  ¥  -  p(|- »)] 


The  following  Information  le  assumed  in  order  to  complete  the 
formulation: 

Propellant  burning  rete^(r): 

t  m  kp-°  ,  w  «  0  }  rdt,  Q  *•  fit tf/w 
o 
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Initial  and  Boundary  Conditions: 

The  gases  are  at  rest  at  the  breech.  The  gases  at  the  base  of 
the  projectile  are  allowed  to  move  with  the  projectile  velocity  because 
of  a  continuity  requirement.  Even  though,  the  continuity  requirement 
is  assumed,  this  is  still  an  unknown  boundary  condition.  However,  a 
relation  can  be  formulated  by  consideration  of  a  dynamic  balance  for 
the  moving  projectile. 

dV 

i.e. ,  m  ^  *  pA-F 

Two  numerical  examples  are  generated,  one  of  which  is  applicable 
to  small  arms  and  the  other  to  an  artillery  weapon.  Further  details 
are  available  elsewhere  [5] .  The  highlights  are  as  follows: 

•  As  the  bullet  moves  ahead  because  of  the  high  pressure  gases,  it 
Issues  rarefaction  waves  which  will  travel  toward  the  breech.  The  re¬ 
flected  waves  from  the  breech  will  interact  with  the  oncoming  rare¬ 
faction  waves  and  thus  form  a  network  similar  to  a  coordinate  grid 
system.  The  path  of  these  waves  may  be  called  characteristic  curves 
(Figure  1  of  reference  5) .  The  applicable  ordinary  differential  equa¬ 
tions  along  these  curves  can  be  solved  by  finite  difference  techniques. 

•  The  density  distribution  it  various  times  is  shown  in  Figure  5  of 
reference  5.  The  pressure  distribution  is  stellar  to  the  density  dis¬ 
tribution.  The  uniform  density  between  bullet  end  breach  is  difficult 
to  justify. 

•  The  velocity  distribution  is  shown  in  Figure  7  of  reference  5* 

For  the  first  portion  of  bullet  travel  (approximately  20  per  cent), 
the  linear  velocity  distribution  is  not  true. 

•  In  general,  the  assumption  of  a  linear  velocity  gradient  le  better 
than  the  assumption  of  uniform  density.  Kotc  that  the  uniform  density 
assumption  tepllee  s  linear  velocity  gradient  due  to  the  equation  of 
continuity, 

•  The  velocity  of  the  gee  at  any  location  increases  with  time  and 
decreases  significantly  even  before  the  bullet  leaves  the  musale. 

•  bines  this  solution  procedure  ie  not  significantly  more  difficult 
than  the  current  practice  **f  living  only  ordinary  differential  equa¬ 
tions  by  consideration  of  apace  mean  averages,  it  la  recommended  to 
consider  at  learnt  the  model  presented  berm  in  say  interior  ballistics 
studies. 
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III  BOUNDARY  LAYERS 

The  separation  of  gas  flow  Into  core  flow  and  boundary  layer 
flow  and  procedures  to  solv Independently  can  be  justified  for  thin 
boundary  layers,  i.e.,  for  low  viscosity  and  large  Reynolds  number 
flows.  Otherwise,  one  may  be  able  to  solve  either  simultaneously  or 
oue  after  the  other  by  some  iteration  scheme.  Under  any  circumstances, 
one  should  have  p  capability  of  handling  laminar  and  turbulent  cases 
of  unsteady,  viscous,  coo}  ressible  flow  with  arbitrary  pressure  gradi¬ 
ents  on  a  bore  surface. 

The  rate  of  heat  transfer  from  the  hot  propellant  gases  to  the 
barrel  is  controlled  by  the  development  of  the  boundary  layers.  The 
flow  in  the  gun  barrel  boundary  layers  could  be  laminar,  transitional, 
or  turbulent  in  nature.  The  type  of  boundary  layer  at  a  particular 
cross  section  at  any  instant  need  not  be  the  same  as  at  another  instant 
Since  the  flow  has  to  start  from  rest  and  must  also  satisfy  zero  boun¬ 
dary  layer  thickness  at  the  bullet  base  hecause  of  the  scraping  action 
of  bullet,  laminar  flow  always  exists  in  some  parte  of  the  gun  barrel 
boundary  layers. 

Flow  in  a  laminar  boundary  layer  will  eventually  become  un¬ 
stable  as  the  Reynolds  number  is  increased.  The  boundary  layer  thick- 
ness,  skin  friction,  and  heat  tranafct  increases  many  fold  for  turbu¬ 
lent  flow  over  laminae  flow.  The  mechanics  of  turbulence  or  more 
common  eddy  viscosity  is  the  dominating  mechanism  for  such  Increases. 
The  boundary  layer  flow  can  be  turbulent  somewhere  in  the  middle  of 
the  flow  between  the  breech  and  the  bullet  base.  transitional  regime 
should  exist  between  the  laminar  and  turbulent  regimes,  however,  be¬ 
cause  of  limited  knowledge  about  transitional  regimes,  the  flow  will 
be  assumed  to  change  suddenly  from  laminar  to  turbulent  at  a  time  and 
a  place  determined  by  a  veil-known  laminar-turbulent  transition  cri¬ 
teria.  Therefore,  the  unsteady  boundary  layer  analysis  is  needed  for 
laminar  and  turbulent  boundary  layers. 

The  state  of  art  in  unsteady  boundary  layers  {6)  and  tur¬ 
bulent  models  f?,6)  is  Quite  limited.  The  literature  cannot  be  re¬ 
viewed  here  due  to  lack  of  space.  However,  YalsaencbiH  {9,10)  re¬ 
viewed  it  in  some  detail  in  the ■ pee t*  It  is  concluded  that  no  other 
investigator  has  solved  an  unsteady  compressible  laminar,  and  turbulent 
boundary  layer  problems  with  arbitrary  pressure  gradients  and  arbitrary 
free  stream  conditions.  The  governing  boundary  layer  equations  with 
eddy  viscosity,  c*,  formulation  in  dimensionless  form  are  as  follows: 

CooctoU,:  If  ♦  £  to.)  ♦  £  (p.)  -  0 
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The  goveruiog  equations  are  a  system  of  nonlinear  parabolic 
partial  differential  equation*  with  three  independent  variable#.  So 
general  solutions  exist  which  allow  arbitrary  initial  and  boundary 
condition#.  Various  solution  procedures  are  attempted  because  of  com¬ 
plexity  and  i«port*nc*.  One  approach  involves  the  use  of  method  of 
weighted  residual#  (method  of  moments)  and  the  method  of  characteristics 
in  order  to  reduce  only  continuity  and  momentum  aquation#  to  ordinary 
differential  equations  and  numerical  integration  in  the  end  (if.  Thus, 
the  u«e  of  existing  analogy  between  moueatu*  end  energy  transfer  for 
steady  flown  ia  required  to  predict  the  convective  heat  transfer,  to 
check  the  validity  of  analogy  for  unsteady  cases ,  another  approach 

•Referenced  quantities 
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involves  the  use  of  some  techniques  on  continuity,  momentum  and  energy 
equations,  but  with  some  more  assumptions.  If  the  analogy,  such  as 
Colburn,  is  valid  for  unsteady  boundary  layers  with  arbitrary  pressure 
gradients,  the  convective  heat  transfer  coefficients  computed  by  two 
approaches  for  any  example  should  be  almost  the  same.  However,  these 
are  not.  The  results  based  on  the  analogy  underestimate  as  much  as  100 
per  cent  for  small  times,  however,  the  same  analogy  overestimate  the 
convective  heat  transfer  coefficients  up  to  50  per  cent  for  large  times. 
It  is  surprising  that  the  Colburn  analogy  predict  almost  identical  re¬ 
sults  as  obtained  without  the  use  of  an  analogy  at  one  of  the  times 
considered.  This  time  is  about  one-half  of  the  projectile  travel  time 
inside  the  tube. 

Since  more  assumptions  are  involved  in  the  above  approaches 
than  desirable,  other  solution  procedures  are  attempted.  The  trans¬ 
verse  coordinate  is  modified  to  absorb  the  compressibility  effect.  The 
stream  function  is  introduced  to  satisfy  the  continuity  and  thus  to 
eliminate  one  of  the  dependent  variables.  In  the  analytical  ap¬ 
proach  [11],  the  method  of  weighted  residuals  is  used  to  reduce  by  one 
the  number  of  independent  variables.  No  way  presently  seems  to  be 
available  to  select  the  functions  systematically  for  all  problems. 
However,  the  functional  in  solution  form  for  the  dependent  variables 
is  chosen  on  the  basis  of  the  asymptotic  solution  of  steady  differen¬ 
tial  equations  for  large  values  of  the  space-like  coordinate.  The 
error  functions  consequently  occur  for  forced  convective  flows.  The 
weighting  function  (Galerkin)  is  chosen  on  the  basis  of  least  error 
for  wall  shear.  All  integrations  across  the  boundary  layer  are  per¬ 
formed  analytically.  The  Method  of  Lines  is  used  to  eliminate  another 
independent  variable.  Finally,  numerical  solutions  are  obtained  by  a 
fourth  order  Runge-Kutta  integration  scheme.  Since  no  real  test  case 
exists  for  such  a  general  solution  procedure,  the  Rayleigh  Blasius  in¬ 
compressible  flow  on  a  plate  was  considered  [9,10].  The  results  are 
in  satisfactory  agreement  with  Hall's  transient  results  and  Blasius 
steady  state  results  with  even  one  term  in  solution  form.  Another 
example,  the  shock  induced  compressible  boundary  layer  problem  is  in¬ 
vestigated  because  of  its  similarity  to  bullet  induced  boundary  layer 
problem.  Again,  the  results  are  satisfactory.' 

The  above  procedure  may  be  tedious  for  turbulent  boundary  layers 
atid  also  for  arbitrary  Initial  and  boundary  conditions.  Therefore,  a 
second-order  accurate  numerical  scheme  is  being  developed  for  the  anal¬ 
ysis  of  laminar  and  turbulent  boundary  layers.  The  scheme  has  many 
desirable  features,  theoretically  and  computationally.  For  example, 
the  Richardson's  extrapolation  scheme  is  adaptable  and  thus  possible 
to  increase  the  accuracy  by  another  two  orders  of  magnitude.  However, 
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the  preliminary  results  indicate  drastic  increase  in  computational 
times  over  the  previously  discussed  techniques  and  also  quite  sensitive 
to  the  input  of  initial  and  boundary  conditions.  But,  this  is  a  valu¬ 
able  tool  to  qualify  analytical  models,  to  study  the  structure  of  un¬ 
steady  boundary  layers,  and  also  to  determine  accurately  the  highly 
unsteady  convective  heat  fluxes. 

IV  HEAT  DIFFUSION 

In  addition  to  various  investigators  all  over  the  world,  Yala- 
manchili  and  Chu  [12,13,14,15]  investigated  finite  element,  finite 
difference,  and  weighted  residuals  methods  for  transient  two-dimension¬ 
al  heat  diffusion  problems.  The  governing  equations,  initial  and 
boundary  conditions,  stability  and  oscillation  characteristics,  and 
applications  are  available  in  those  references. 

V  FREE  CONVECTION  AND  RADIATION 

The  outer  surface  temperature  of  the  gun  tube  may  change  from 
-60°F  to  1500°F  because  oi  the  high  rate  of  fire.  Since  wind  veloci¬ 
ties  help  to  cool  the  gun  tube  faster  than  free  convection  and  radia- 
tic  „  thermal  design  of  a  gun  tube  should  not  Include  wind  velocities. 
Tb:'  ,  utation  of  the  product  of  the  Grashof  and  the  Pr&ndtl  numbers 
for  typical  gun  barrel  conditions  (106)  indicate  that  the  flow  is  in 
the  laminar  regime  (10 3  to  109).  The  contributions  of  free  convection 
and  radiation  are  almost  equal  where  the  outer  surface  temperature  of 
the  gun  tube  is  400°F.  The  free  convective  contribution  dominates  over 
radiation  for  temperatures  much  less  than  400°F  and  vice  versa.  The 
governing  equations  of  unsteady  free  convection  and  radiation  are  a 
system  of  nonlinear  partial  differential  equations  of  parabolic  type 
with  three  independent  variables  and  with  variable  linear  and  non¬ 
linear  (radiation)  boundary  conditions.  Yalamanchlli  [16]  discussed 
the  solution  procedure  by  explicit  finite  differences  and  also  results. 
The  conclusions  are  as  follows:  The  results  are  stable  even  with  40 
per  cent  larger  time  step  sizes  than  obtained  by  stability  criteria. 

The  thermal  boundary  layer  in  the  upper  half  of  the  cylinder  is  much 
thicker  than  in  the  lower  half  of  the  cylinder.  The  thermal  boundary 
layer  thickness  increases  not  only  with  tine  but  also  with  position 
for  large  times.  The  velocity  boundary  layer  thickness  also  grows 
with  time  but  not  as  fast  as  the  growth  of  the  thermal  boundary  layer 
thickness.  The  dimensionless  heat  transfer  is  fairly  uniform  for 
small  times  and  changes  up  to  50  per  cent  over  its  maximum  at  that  time. 

COMBINED  ANALYSIS 

The  solutions  obtained  for  the  individual  fictitious  problems 
do  not  represent  the  solutions  for  the  real  gun  tube  problems  because 
of  continuous  change  in  boundary  conditions.  Compatible  boundary  con¬ 
ditions  must  be  introduced  at  Interfaces  between  the  problems.  It  is 
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convenient  to  generate  ahead  the  thermochemical  properties  for  various 
pressures  for  any  chosen  propellant  and  later  input  in  the  form  of 
tables  for  problems  II  and  III.  To  solve  core  flow  problem,  one  may 
require  the  location  of  boundaries,  in  particular,  the  interface  be¬ 
tween  problems  II  and  III.  Since  the  location  of  this  interface  is 
unknown  and  also  the  thickness  of  boundary  layers,  in  general,  is 
small,  one  can  reasonably  assume  that  the  interface  is  the  bore  surface 
for  problem  II.  However,  one  should  include  the  information  at  the 
interface  such  as  the  heat  loss  to  the  gun  tube  and  skin  friction. 

These  quantities  will  depend  on  the  development  of  boundary  layers  and 
also  the  temperature  of  the  bore  surface.  Sometimes,  the  assumption 
that  the  outer  edge  of  the  boundary  layer  extends  to  the  centerline  of 
the  tube  because  of  negligible  viscous  and  hast  conduction  effects  in 
core  flow  is  convenient. 

Even  though  the  boundaries  are  well  defined  for  a  chosen  gun 
for  heat  diffusion,  the  information  at  the  boundaries  such  as  heat-in 
due  to  forced  convection  and  heat-out  due  to  free  convection  and 
radiation  are  lacking  because  these  in  turn  are  dependent  upon  pro¬ 
blem  IV.  It  seems  that  there  is  strong  coupling  between  boundary  layer 
(forced  convection)  problem  and  heat  diffusion  problem  because  of  rapid 
change  in  conditions  at  the  Interface.  The  combined  unsteady  forced 
convective  and  heat  conductive  problem  will  be  solved  once  analytical 
capability  in  unsteady  boundary  layers  is  well  established.  The 
classical  heat  diffusion  equation  possesses  several  pathological 
anomalies,  the  most  important  one  being  an  infinite  speed  of  heat  wave 
or  thermal  shock  propagation.  The  use  of  the  Fourier  model  is  definite¬ 
ly  questionable  for  heat  pulse  problems,  plasma  torch,  etc.,  where  high 
heat  fluxes  may  be  present.  Therefore,  a  non-Fourier  model  Is  being 
investigated  because  of  the  possibility  that  the  bore  surface  tempera¬ 
ture  and  its  gradient  could  be  even  higher  than  the  classical  results. 
The  heat  transfer  rates  due  to  frictional  heating  between  rotating 
band  and  rifling  at  the  rubbing  contact  are  very  high.  Probably,  these 
could  be  much  higher  than  the  propellant  gas  convective  heat  transfer 
rates.  Similarly,  the  forced  convective  heat  transfer  will  be  sub¬ 
stantially  higher  for  hypervelocity  systems  because  of  higher  gas  velo¬ 
cities  and  also  large  temperature  differences  between  the  gas  and  the 
tube.  The  effect  of  non-Fourier  effects  could  be  significant  at  least 
in  such  circumstances.  Finally,  it  should  not  be  forgotten  about  the 
possibility  of  interaction  between  core  flow  and  boundary  layer  flow 
for  hypervelocity  systems  such  as  In  hypersonic  flow  [17].  It  is  also 
Important  to  recognize  the  existence  of  shock  ahead  of  the  projectile 
for  hypervelocity  systems  and  consider  its  effects,  such  as  resistance 
induced  on  projectiles  (core  flow)  and  also  shock  induced  compressible 
boundary  layers.  All  the  models  proposed  are  quite  rigorous,  but  fail 
to  forecast  why  the  overall  objective  can  not  be  fulfilled. 
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